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Abstract 

We develop the necessary tools for computing fluctuations around a mean-field background 
in the context of SU(N) lattice gauge theories in five dimensions. In particular, expressions for 
the scalar observable and the Wilson Loop are given. As an application, using these observables 
we compute a certain quantity k§ that can be viewed as Coulomb's constant in five dimensions. 
We show that this quantity becomes independent of iV in the large N limit. Furthermore, the 
numerical value of we find for SU(oo) deviates by 17% from its value predicted by holography. 



1 Introduction 



In four dimensions and in infinite volume SU(N) gauge interactions are renormalizable and 
always confining. Five-dimensional (5d) SU(N) gauge theories are non-renormalizable and 
trivial at the perturbative limit. Their weak coupling phase is Coulomb and they have a first 
order phase transition at a critical value of the coupling. Beyond the phase transition there 
is a confined phase. The appropriate analytical method to describe these theories near the 
phase transition from the side of the Coulomb phase is the Mean-Field (MF) expansion [1]. 
The motivation to consider these theories, beyond the possibility of the existence of a physical 
fifth dimension, is that the MF method becomes more accurate as the dimension of space- 
time increases. This gives us the opportunity to develop a trustable analytical probe away 
from the perturbative point. Then one can look at regimes where the system is dimensionally 
reduced and build effectively four-dimensional (4d) models. This would not be possible directly 
in four dimensions as the MF background vanishes identically in the confined phase. As we 
argue here however, there are interesting theoretical issues that can be discussed already in 
five dimensions, in particular related to the large A limit of these theories and their possible 
holographic description. 

In [2] the five-dimensional SU(2) MF formalism was developed, formulated on a periodic 
and anisotropic lattice. We will be using this work's notation and results extensively. In a first 
application [3] it was shown that there is a regime on the phase diagram where the system reduces 
dimensionally by localizing the gauge interactions on four-dimensional hyperplanes. Interactions 
on the hyperplanes are confining and the string tension as well as some of its corrections, such 
as the Luscher term, were computed. The localization of the gauge interactions in this regime 
was consequently supported also by Monte Carlo simulations in [4]. In the meantime there has 
been a steadily increasing activity in lattice simulations of five dimensional gauge theories [5]. 
In [6] the SU(2) MF model with orbifold boundary conditions is constructed, the goal being a 
non-perturbative description of Gauge-Higgs Unification. These (semi) analytical computations 
complement earlier lattice orbifold Monte Carlo simulations [7]. 

Another application of the MF formalism, not relying on dimensional reduction and one 
which will be of our main interest here, was presented in [8] where a certain quantity, an analogue 
of Coulomb's constant in five dimensions was computed. It is derived from the five-dimensional 
Coulomb static potential Vs(r) = const. — C2/r 2 , where C2 is the static charge, as follows: 

h = c 2 ^, (1) 

where C2 = C2/a, with a the lattice spacing and (3 the lattice coupling. The result, obtained 
from SU{2) gauge theory on periodic and isotropic lattices, was compared to a holographic 
computation of the same quantity [9]. There it was argued that this quantity, at least in the 
large A limit, is in fact A-independent and this could justify a comparison of the A = 2 lattice 
result with the large A gravity result. The two calculations of k$ showed a numerical agreement 
of 2% [8], suggesting that indeed, even if k*, is A-dependent, its dependence on A is probably 
rather weak. 
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In this work we extend the SU(2) MF formalism of [2, 3] to SU(N). In [1] one can find 
expressions for the SU(N) MF propagator on the isotropic lattice in the axial gauge. Here we 
first extend the results of [1] by computing the propagator in a covariant gauge on the anisotropic 
lattice. Moreover, we compute the Wilson Loop and the scalar observables, that have not been 
computed before for SU(N). As a first application, these quantities allow us to compute k$ for 
large N and verify its weak iV-dependence suggested by holography. 



2 Mean-field formalism 

The general formalism of the MF expansion was presented in detail in [2, 6] and we will not 
repeat it here. For details one can consult these references together with [1]. One is interested 
in physical observables, schematically written as 

(0} = ^ J BUO[U}e' Sw ^ . (2) 

Here Sw [U] is the usual Wilson plaquette action with anisotropic couplings /?4 = f3/j along four- 
dimensional hyperplanes and /?5 = /?7 along the fiith dimension. In this notation /3 is referred 
to as the lattice coupling and 7 as the anisotropy parameter. The isotropic case is obtained for 
7 = 1. To first order in the MF expansion the expectation value of an observable takes the form 



with 



_K- 1 , (3) 
v J 



K = _ K {hh) 1 + K (w) + K (gf) (4) 

the lattice propagator. K^ hh ^ is the second derivative of the MF effective action with respect 
to auxiliary degrees of freedom H and K is the second derivative of the Wilson plaquette 
action with respect to the link variables V. is the gauge fixing term. O is a gauge 

invariant operator and its second derivative with respect to the .ff-variables contracted against 
the propagator and evaluated in the MF background defines an expectation value. The connected 
version of the latter is defined as O c (t) = 0{t^ + t)0{to) and from it the correlator 

C(t) =< O c {t) >-< O(t + t) >< O(t ) > (5) 

is formed. The ground state mass associated with the operator O, to first order in the fluctua- 
tions, is then extracted from 

C(t) 

2.1 Some SU(N) integrals 

Here we review some facts related to SU(N) integrals [1]. Let f{U) be a function of the SU(N) 
link variables U such that f(U) = /(VUV- 1 ), V G SU(N). Also let X r(U) be the character 
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associated with the irreducible representation r of SU (N) that the group element U is expressed 
in. Then, 

f(U) =Y,frXr(U), f r = [ DU f{U) X *{U) (7) 

r J 

A special case is the function 

e h l^)=J2crXr(U). (8) 
r 

We express the irrep r in terms of a set of positive integers as [10] 

r -> {A«} = {A« > A< r) > A« > ■ ■ • > A« = 0} . (9) 

We will need the A^ r ^ only for the symmetric and anti-symmetric representations. They are: 

A< 2 > = (2,0,--- ,0,0) 
A (1 ' 1} = (1,1,0,--- ,0,0) (10) 

respectively. Then, 

e > etr([7) = E X{x(r)}[u) g det U (r)j _ j+i+n (M } , 

|^(r)| n=-oo ^ \ / ) 

with i,j = 1, ■ ■ ■ , N and I q the Bessel function of order q. Using that 

J DU X{X (r )} (U) = 6 {x(r)}m (12) 



we have 



fAfVn ra=— no \ / 



{AM} 
+oo 

= £ det j J_ i+i+B ( ^ ) \. (13) 



+°° r / 1 



An example is N = 2 where 



/ Due**** 1 * = ^det ( In . (14) 

J n \In+l In J 

In a matrix representation we can write a link in the irrep {A^ r )} as the matrix 

U(n',M')^D { Jp(n',M'). (15) 
The character of {A^} can then be written as 

X { xtr)}(U) = D^ r) y . (16) 
We note the useful orthogonality relation: 

/ DUD%^D& r >> = S {xlr)ux{rY} S^^- • (17) 
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2.2 K ihh ^ 



Let U G SU(N) the link (with explicit indices U a p(n,M)) and H an N x N complex matrix. 
The indices a, /3 = 1, . . . , N are the gauge indices, n is the location of the link in the lattice and 
M is its direction. Define 

' DUe^ Kctr{ - UH) 



((H) = 

u(H) = -logC(tf). (18) 

We are after derivatives of u{H) with respect to H a >p/(n' , M'), evaluated in the mean-field 
background 

(19) 



V = VI, 



H = HI 



with 1 the N x N unit matrix. We define the basic integral 

+oo 



Co = / DUe****™ = £ det {/_ J+l+n (£) } 



(20) 



where in a given parametrization V = vq and H = Hq. We start from the first derivative of 
u{H) 

du{H) i ac(-ff) 



<V/3' - dHp, a ,{n>, AT) lH=H1 C dHpwiri, M>) lH=H1 



1 1 



DU 



d 



dHp, a ,{n',M') 



1 

N 
1 1 



Re tr([/F) 



^Retr([/fT)| 



\H=H1 



(21) 



where 



Ci(a , ,/3 , ,n',M') = J DU U a ,p,(n' , M') e^^( u ^ 5 a ,^ 

Ci = / dU(TrU)e^ 
JsU(N) 



'-ReTrU 



(22) 



In the above expression we have dropped the argument (n',M') from Ci«'/3' because the back- 
ground is uniform. In the anisotropic background the effect of the anisotropy is encoded in 
the vev H (it will be ho along the /x-directions and ho5 along the fifth dimension). The second 
derivative of u(H) is 

d d 



u a',l3' ] a",l3" = 



MH)\h- 



HI 



dHp W (n>,M>)dHp //a „(n»,M") 

[Ci(a', /?', n', M')Ci(a", /?", n", M") - CoC 2 (a', /?', n', M'; a", /3", n", M")] 



with 



( 2 (a,(3',n',M';a",(3",n",M") 
d 



(23) 



/ 



DU 



_dHp, a/ {n',M>) 



Re tr([/tf) 







dHf3„ a „{n",M 



TTTUr^e tr(UH) 



±R C tv(UH) 



H=H1 



(24) 
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that is, 



( 2 (a',(3',n',M';a",(3",n",M") = J DU U a ip(n' , M 1 ) U a »pi>(n", M") e 7v Retr ( f/ ) . (25) 

Next, we compute the integrals in Eq. (22) and Eq. (25) respectively. We will make repeated 
use of the notations: 



D n ,N{z) = det [In+j-i ( z ))l<i,j<N ' D n = D u,n(z 



D 



, dD ntN (z) 



dz 



D 



„ _ d 2 D ntN (z) 



dz 2 



.ho 

' N 



A kind of generating function for SU(N) has been presented in [11]. It reads: 



G 



SU{N)\ {C;d) 



Isi 



^jj e cTrU+dTrU^ 

'SU(N) 

We note the following consequences: 



E 



det 



l n+j- 



(2\fcdj 



l<i,j<N 



G 



SU(N)\(hQ hp 



dG 



SU(N) 



dc 



dG 



SU(N) 



d 2 G 



dd 

SU(N) 



dc 2 



( ho_ h o \ 

\2N ' 2N J 

( H. h o \ 

\ 2N ' 2N J 

I h^_ hp \ 
\2N'2N J 

( ho_ h () \ 

\2N ' 2N J 



d 2 G 



SU(N) 



dd 2 



\2N ' 2N J 



[ dUe^ ReTrU = V det 
J SU(N) n= _ oc 

r h +°° 

/ dU(TrU)e^ ReTrU = V 
JSU(N) 

/ dU(TrU^ ReTrU = Y 
JSU(N) n ±^ c 

[ dU{TrU) 2 e^ 
JSU(N) 

+oo 

E 

n=— c 

/ dU{TrU ] ) 2 e^ 

JSU(N) 



ln+j-i 



l<i,j<N 



V » h Ur 



nN 2 ^ 



h 



ReTrU 



D»--(l + 2nN)iy n+ ^ 



ReTrU 



2nN 3 fnN 



+ 1 A, 



'SU(N) 
+oo 



= E 



n h [ 



, . 2nN 3 fnN \ 
1 - 2nN)D' n + - l j D r , 



(26) 
(27) 

(28) 



d 2 G s , 



U(N) 



dcdd 



( h. h o \ 

\ 2N ' 2N J 



[ dU(TrU)(TrU^)e^ 
JsU(N) 



ReTrU 



E 



N , n 2 N 4 
n" j n' n 

n+ ho n h 2 Un 



(29) 



Since D n is even under n —> — n, all terms that are odd in n vanish in the sums. For the integral 
Ci(a', /3', n', M') we have 

1 



Ci(a',(3',n',M') 
Ci 



N 



SU(N) 



dU[TrU(n',M')]eW ReTrU = %f Ci, 



AT 



+oo 

E 



(30) 
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Thus, in this notation, 

oo oo 

co= Yl Ci= E D 'n- ( 31 ) 

n=— oo n=— oo 

For the integral £2 we have to distinguish three different cases. The first is when we take two 
derivatives with respect to H' and H" , the second is when we take one derivative with respect 
to H 1 and one with respect to H"* and the third is when we take derivatives with respect to 
ff't and £T"t. We call the corresponding integrals as (2°, C° + (C^ = (C2 + ) f ) and C 2 f+ (= £2°) 
respectively. We take all iJ's in the fundamental representation. Again, we can drop all space- 
time and directional arguments since all the information about the background is contained in 
the exponent of the integrand in Eq. (25). We have 



i 



V. ,V, M<; a", „», M<<) = / DU V„ K, M') („», M<<) - MP) (32) 



Ci°5 q '/3'<5q;"/3" + C2°5 a ipti5 a iipi. (33) 



We have set the integral equal to the two possible tensor structures, with coefficients to be 
determined. Contracting with S a >p/5 a »pn both sides we obtain 



/ 



DU [ X (l0] 2 e> c tr(l/) = c?°iV 2 + cfN ee J^, (34) 
while contracting with 5 a ip"5 a "pi we obtain 

£[/ X (^)e> ct ^ = c ?°iV + 4°N" = I a p,p a . (35) 



/ 



Thus the coefficients can be expressed in terms of the integrals: 

no _ N J a«;/3/3 ~ tgfcPa no _ jY Ia ^ a ~ ,o R \ 
Cl " N(N*-1) ' 2 " N(N*-1) ' W 

Following [11] we find that 

W = / ^Mt/) + xii(£/)]e> etr ^ 



+00 / h \ + °° 

- E ^ ^ E de«/ 

n=— 00 v 7 n=— 00 



ho 

j-i+n+\f 1) \ ]Y 



(37) 



- E <•* W»?» ( h 4) - E w*p» (v ) ■ < 38 > 



j-i+n+AV ^ \ N J ^ j-i+n+\\ ' V _/V 

rt=— 00 v 7 n=— oo v 

These integrals will be computed numerically. The ++ case does not require extra work, we 
have that cf + = c? and c^~ + = c!] . The 0+ case needs though separate computation. The 
integral to be computed is 

ti 0+ \a'f3';a"(3") = J a ,p,., a „p„ = J DUU a ,p,(n> ',M')^„(n" M") e^ 6 ^ (39) 
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Following the same steps as before we find: 

c^N 2 + 4 + N = J aa . M = f DUtr(UMU^)e^^ u \ 

c°+N + = J a ^ a = I DC/tr([/[/t)e>^) = N( , 



(40) 



which yields: 



0+ Jaa;(S/3 ]\[Ja[3;f3a 



„0+ 



^—^ J DUe^***™ (TrUTrrf - l) , 



(41) 



(42) 



(43) 



We can finally write down K^ hh ^ directly in momentum space, since the Fourier transformation 
is trivial (see [2]). We introduce the index q', q" = 0, +. Then, 



K^ MII {p',a'^', q '-p'',a''^\ q '') 



~ 6p,p "WCo 

= ~ S P'P" WCo 

6p,p "WCo 



CI 9 (a>,P';a",l3") 



N 2 (o 



v ,h 



diag(l,l,l,l,0) 



N 2 (o 



V 05,^05 



diag(0,0,0,0,l) 



2 Sa',P'S a ",P" 



v ,h 



I q'q" x x , i'q"x x ^ z-2 S a',p'S a ",p" 

(q a 'P'0 a "P" + C 2 a '/3"d a "j3>) - Ci 77^7 

^ Co 



^05,^05 



•diag(l,l,l,l,0) 

diag(0,0,0,0,l). 

(44) 



The Euclidean indices M'M" are shown explicitly in (diagonal) matrix form. Before we continue, 
we will rewrite this expression in a more useful basis. To begin, evaluating second derivatives 
in the mean- field background and leaving the gauge index structure aside for the moment, 
corresponds to 



(00): Re[d u d u ] = l/4(d 2 UQ -d 2 UA ) 
(0+): Re[d u d u *] = l/4(d 2 U0 +d 2 UA ) 



(45) 



where the link is parametrized as U = uq+iua- Therefore the Hermitian (H) and anti-Hermitian 
(AH) channels are 



H = (00) + (0+) 
AH = (00) - (0+) 



(46) 



On the other hand, in order to handle the gauge structure, one introduces the projectors on the 
singlet and adjoint representations 

1 



P (S) 
p(A) 



N 



8 a >p"5 a "p" 



8a'P"8 a "P' - -Jj&a'P'ba"P" 



(47) 
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and the new coefficients 

cf = 4 00) + 4° +) 

= c^-c^) (48) 

with a = 1 , 2 and 

A ~ iv 2 Co (c2 +iVci W' "iv 2 Co C2 
^ = -]^(^ H + ^). bAH = Wco C * H (49) 

Now we can write the non-vanishing components of K^ hh ^ simply as 



vo,ho 

Ki h 5 h \p',a',P';p",a",P";z) = 6„„. >U >P™ ;Q „ ^ B* )_ (50) 



where z = H, AH. One can check that for any iV 



Nh Nh xN$ Co(x) ' 

We have defined x = Hq/N and used the background solution given later in Eq. (67). 

For large N, the numerical computation of the determinants involved in the coefficients is 
plagued by instabilities due to large cancellations. It is then useful to consider the large H 
expansion of u(H) [1] 



u(H) = In 



(2tt)( n - 1 )/ 2 



N 2 -2 - N 2 - 1 - iV 2 - 1 

+ InN + H-- InH- . +••• (52) 

2 2 8H K ' 



with H = tr ^^ , which is equal to ho in the mean- field background. Noticing in addition that 
as N increases (3 C and ho also increase we can conclude that already for relatively low values of 
N (N = 5, 6, • • • ) the asymptotic expansion of u(H) is a good approximation. 

In this limit we can use the approximate expressions 

. JV rr 

K 

A AH ~ 0, B AH = _^ (53) 

Nho 

and for later reference 

™ = ^ = 1 -^sT' (54) 

For N = 2, B H and A AH are identically zero and for large N, they become negligible [1]. Notice 
also that the expression for B AH is exact for any N, see Eq. (51). For instance, on the isotropic 
torus we have 

1 1 

When ho — > oo we have vq ~ 1 and we obtain the asymptotic expression B AH = — - 



bAH = -it^=s- ( 55 ) 



Nho 



2.3 



The space-time structure of for SU(N) is similar to the SU(2) case. In q' — q" space it 

is a two by two matrix, just like K^ hh ^ . The hermitian (anti-hermitian) channel of the SU(N) 
model corresponds to the hermitian (anti-hermitian) channel of the SU{2) model. The gauge 
index structures of and of the gauge fixing term can be seen from the second derivative 

of the Wilson plaquette action (suppressing all but the group indices) 



d 2 tr 



u aP u M u\ s ul 



and the gauge fixing term 



dUa'p'dUa'/p" 

9U a i pi dU a » pa 



5 a 'p"S a "i3'N 



respectively. 

We define the bond shifting operators 



( 



7 

Y' c m> - 

ySl/2S0/2 
2/ s 2/2 s 0/2 
2/S3/2S0/2 
\ y5S5/2«0/27 2 ^ 



±s 2 



2/SO/2S1/2 



is 2 



Y' C M' - 
U s 2/2 s l/2 
2/S3/2S1/2 



y5S 5 /2Sl/27 2 ^ 



VS0/2S2/2 
VS1/2S2/2 

El 1 2 

CM' - I s 2/2 

2/S3/2S2/2 
y5S 5 /2S2/27 2 |f 



2/S0/2S3/2 
VS1/2H/2 
2/ s 2/2 s 3/2 

El 1 2 

CM' ~ 

2/5S5/2S3/27 



2 VQ5 
VQ 



where 



2/5 = 27 



2^05 7 



VQ 



and 



A„ = -(iV^ 2 )- 

/ I]' CM' 

2Ci/2Co/2 
2c 2 /2Co/2 
2C3/2C0/2 

V 2c 5/2 c 0/ 27 2 t 



(56) 



(57) 



y5S0/2«5/27 2 ^ 
y5Sl/2«5/27 2 ^ 
y5S 2 /2S5/27 2 lf 
y5S3/2«5/27 2 U 



7' 



CM' 



1,2 



(58) 



(59) 



2c /2Cl/2 
I]' C M' 
2c 2 /2Cl/2 
2c 3 / 2 Ci/2 

2c 5 / 2 Ci/27 



2c / 2 c 2 /2 
2C1/2C2/2 

Y' C M> 

2c 3 /2C 2 /2 



2C0/2C3/2 
2C1/2C3/2 
2C2/2C3/2 
Yl C M> 



2 VQ5 



2c 5/2 C2/27 2 |f 



2c /2C 5 /27 2 |f 
2c 1/2 c 5/2 7 2 |f 
2c2 /2 c 5 /27 2 |f 
2c 3 / 2 C5/27 2 ^ 

2c 5 /2C3/27 2 |f 



1 2 Y'CM> J 



(60) 



We use the notation s / 2 = sinp' /2 etc. and C5 = 7 2 ^cos(p' 5 ). These expressions appear in 
[1] in the axial gauge and for isotropic lattices. Here we have generalized them to anisotropic 
lattices and we have fixed a covariant gauge parametrized by £ as in [2]. We have checked 
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numerically that our results do not depend on the value of £. Using 8 a ipn8 a iipi = P^ + P^ 
we can express as 

K M>M"(p'i Q '' Z 3 '^ a "> Z 3 "; = S P'p"( P ^l;a"l3" © P i^';a"l3") A HM'M" 

K M>M"(P'' a '> Z 3 '^"' Z 3 "; = V( P ^j'; a "/3" © Pa^- a "/3") A AHM'M" ■ (61) 

The 5[7 (iV) propagator is then 

K M 'M" = ~ K WM" + K M^M" ■ ( 62 ) 

The notation here and from now on is that K^, M „ represents the M'M"'th element of K —1 . 
The inverse of K has four sectors, labeled by z = H, AH and w = S,A. Schematically we 
represent these sectors as K~ 1 (z,w) and express them in this basis as 

K~\z, w) = K~\H, S)PW + K-\H, A)PW + K~\AH, S)P {S) + K~\AH, A)pW . (63) 
2.4 The phase diagram 

In the anisotropic theory, there are two mean values for the links, vo and ^05 determined by the 
extremization of (d = 5) 

Sefi[ ^ H] = -fa {d - ~ 2) 4~ - l)tt§4s + (d- l)u(h )+u(h 05 ) + (d- l)h v + W05 . 

(64) 

Here TV = L Z N$ is the number of spatial lattice points with N5 the number of points in the fifth 
dimension. Eq. (64) yields conditions identical in form to the 577(2) case [2] 

— — 6/3 

v = -u(h )' , h = — vl + 2/3-fv vl 5 , 

7 

^05 = -u(ho 5 )' , /105 = 8/37^0^05 , (65) 
where u', according to Eq. (21) is (here we differentiate form N = 2) 

u' = tv{u a ,p,} = -^N( 1 , (66) 

that is 

"° = Jv&w ■ < 67 > 

For N = 2 this reduces to the known result vq = h(2x) / I\(2x) . 

One can immediately see that when all momenta vanish, in the (AH, A) sector we have (for 
simplicity we take 7 = 1), using Eq. (51): 



K{AH,A) = -(^+AN(fi%yi 



(68) 



which vanishes by Eq. (65). These are the five expected torons, which persist even when 7 / 1. 
The various phases on the phase diagram are defined as follows: 
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• Confined phase: Do = 0,vo5 = 

• Coulomb phase: vq / 0, U05 / 

• Layered phase: vq 7^ 0,tJo5 = 

We will perform our computations in the Coulomb phase where the background is non-vanishing 
everywhere. In addition we will stay near the phase transition where we expect that cut-off 
effects are suppressed. The critical value of the coupling /3 C that signals the end of the Coulomb 
phase moves towards to larger values as N increases. Apart from that, the phase diagram is 
qualitatively similar to the SU{2) phase diagram as described in [2]. 

3 Observables 

To first order in the fluctuations we will be computing 

{0) = 0[V] + \tr{Y J E (S)^-^^)^ 1 )^^)}- ( 69 ) 

w=S,A z=H,AH V ° 

We will often call the above second derivative of the observable simply as 2 which also can be 
expressed as 

2 (z) = 2 {z, S)PW + 2 (z, A)PW . (70) 

Then we can expand as 

(O) = 0[V} + ^tr{0 2 (z,w)K- 1 (z,w)} 

z,w 

= o[V} + 

Ui [0 2 (H, S)K-\H, S) + 2 {AH, S)R- 1 (AH, S)] tr{l} s 

^tr [0 2 (H, A)K~\H, A) + 2 (AH, A)K~ 1 (AH, A)] tr{l} A 
= 0[V} + 

[0 2 (H, S)K-\H, S) + 2 (AH, S)R- 1 (AH, S)] N 
hi [0 2 (H, A)K-\H, A) + 2 (AH, A)K~ 1 (AH, A)} N(N 2 - 1) . (71) 

3.1 The scalar mass 

The Euclidean structure is identical to that of the SU (2) model. We recall the SU (2) result (for 
SU(2) the gauge index takes the values a' = 0, 1,2,3 where oi = represents the Hermitian- 
Singlet channel and a' = 1,2,3 the anti-Hermitian- Adjoint channel) where only the Hemitian 
channel, free of zero modes, is present [2] 

C S (t) = |(Po) 2 £ cos(p' t) £ \AW(iQ\*K£ ((pf,, 0,p' 5 ), a' = 0; (p' , 0,p' 5 ), a' = 0) . (72) 
Po P5 
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Recall that for SU(2) for any observable there are at most the two sectors, (H, S) and (AH, A) 
contributing. 

For SU(N) the contribution of the adjoint channel vanishes identically for basically the same 
reason as for SU(2), that is because taking only one derivative of the observable with respect 
to a fluctuation along the Lie algebra and evaluating it in the background gives O2 (H, A) = 
2 (AH, A) = in Eq. (71). Also, since the double derivative on the Polyakov loop that repre- 
sents the scalar has the gauge index structure 

6a>tl»6a»,fi> = P (S) + P (A) (73) 

and the schematic structure 

C s (t) ~ NQ 2 (H, S)K-\H, S) + N0 2 (AH, S)K~\AH, S) , (74) 

the final expression for the correlator becomes 

NN 2 ( P ) 2 r 
Cs{t) = a/ J E cos (Po*)V fc ,o^,o K£ {p'-p'-H,S)+K£ (p';p';AH,S) 



Po 



(75) 



where k = 1,2,3. Notice that for SU(2) the coefficient A AH vanishes identically, so only the 
(H, S) sector contributes, consistently with our comment above. We have defined 

A^(p) = — Y e^+Va) = 7v 5 fe , (76) 

where r = labels the first link along the fifth dimension, r = 1 the second, etc. Po is the 
mean-field length of the Polyakov loop representing the scalar observable. It is easy to see that 
the toron p^ = adds a zero to the scalar correlator since it contributes a time independent 
constant. Given that the Hermitian channel of the propagator does not contain any zero modes, 
there is no toron. The plateau in the time decay of this observable yields the mass of the mass 
of the scalar observable 047715 in lattice units. 

The vector mass is a purely finite volume effect (in particular it does not depend on N) so 
we can use directly the SU(2) relation [2, 3] 

47T 

a^vriv = — . (77) 
1j 



3.2 The Wilson Loop 

Here we compute the Wilson Loop (WL) along four-dimensional hyperplanes. Again, the space- 
time structure is the same as for SU(2), see [2]. 

For an observable represented by a single loop it makes a difference whether the two links 
that are removed (by the two derivatives) are pointing in the same or opposite directions. In the 
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former case there is a relative minus sign with respect to the latter. A case where this remark 
applies is the Wilson Loop. The schematic structure in this case is 

exchange : (00) + (0+) 
tadpole : (00) - (0+) (78) 

corresponding to the two diagrams that contribute to leading non-trivial order, an exchange 
between the two spatial legs and the tadpole on a given spatial leg. For general N, in principle, 
all four sectors contribute. The final result is: 

° k=i P > 

| [(^cos^r) + l)K^{p';p';H,S) + (N 2 - cos(^r) - l)K^(p'-p'-AH,A) 

"(I cos{p ' kr) _ l) K ^(p';p'; AH, S) + (N 2 - 1)(1 cos^r) + l)K 00 \ P ';p'; H, A)] j 

(79) 

where we have dropped the irrelevant additive constant originating from the zeroth order con- 
tribution. 

For SU(2) only the first two sectors contribute. In particular, the two terms in the last line 
vanish because A AH = B H = 0. Finally the toron in the propagator in the {AH, A) sector is 
cancelled by the zero in the observable in the same sector since the factor cos(p' fc r) — 1 vanishes 
for p = 0. Notice that for SU(2) the coefficient A AH is not being used by neither the scalar nor 
by the WL. 



4 A simple application: Coulomb's constant in five dimensions 

We present a simple application of the formalism we developed, on the isotropic lattice. As 
mentioned in the Introduction, the dimensionless quantity k§ defined in Eq. (1) is computable 
both in lattice gauge theory and in gravity. An agreement between the two computations would 
constitute evidence for the validity of the holographic conjecture [12]. We remind that the 
holographic computation of k$ defined in Eq. (1) yields [9, 8] (expressed in lattice parameters) 
the number 



5(2/3, 1/2) 



0.0649 (80) 



3vr 2 / 3 

with B(x,y) the Euler Beta function, for any N for which the computation is valid. This is 
expected to be the case in the large N limit, therefore any agreement with the lattice computation 
of the same quantity is expected to occur at least in this limit. 

From the point of view of the lattice, we define trajectories on the phase diagram that 
approach the bulk first order phase transition, with 

q = ^ (81) 
a^ms 
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kept constant [3]. The quantity k§ = C2P/N 2 on the isotropic lattice depends on (3 and L. Using 
these trajectories, it is however straightforward to obtain its infinite L value via extrapolation. 
Here we choose q = 2 even though the infinite L extrapolations of k$ are expected to be essentially 
(/-independent. 

In figure 1 we compare for SU(5) the results obtained from the asymptotic expansions 
Eq. (53) for the coefficients in K^ hh ^ together with the asymptotic form for the background 
Eq. (54) against the ones obtained by the sum over determinants Eq. (49). We expect that the 
former (latter) would be more reliable for large (small) values of N and indeed for N = 5 the 
discrepancy between the two is already quite small. According to this result, for iV < 5 we will 
be computing the coefficients using the sum over determinants method and for iV > 5 using 
their asymptotic expressions. 

In figure 2 we show the values obtained for £5 for increasingly large values of N. For 
N = 2, 3, 4, 5 we use the sum over determinants expressions and for iV > 5 we use their asymp- 
totic forms. We also show quadratic fits to the data. We observe that the intercept increases 
slightly with N. The largest value depicted is N = 100. Data for larger N than that become 
indistinguishable from the latter with naked eye. It is clear that the curves saturate for large 
N. This demonstrates that £5 becomes ^-independent as iV — > 00. 

In figure 3 we depict the N dependence of the intercept for N > 10 along with a quadratic 
fit. It cuts the vertical axis at k$ = 0.0757, presumably the result for SU(co). This should be 
compared against the holographic result 0.0649, resulting into a 17% disagreement. Evidently 
the N = 00 value of k$ is farther away compared to its N = 2 value which is just 2% away from 
the holographic result. This "reverse" trend could be related to the order of the computations 
on both the gauge theory and gravity sides. The essential fact to keep in mind is that the 
mean-field prediction for k% in the large N limit stabilizes to a value not too far away from the 
holographic prediction. 

5 Discussion 

In order that the comparison between the Mean-Field and Holographic computations makes 
sense, both should be reasonably good descriptions of the 5d SU(N) gauge theory at intermediate 
to strong coupling. Assuming that from the holographic point of view this is a conjecture and 
provided that the Mean-Field is a good analytical tool at strong coupling, any agreement with 
the gravity side gives evidence for the validity of the holographic conjecture. Reversely, if 
we assume that the holographic duality holds, then any observable derived from the classical 
gravitational system should be in principle verifiable from the gauge theory side, provided that it 
can be computed reliably at the non-perturbative level. No matter what point of view is taken, 
a necessary condition is that the truncated Mean-Field expansion we have employed here is a 
good approximation to the non-perturbative system. We have mentioned in the Introduction 
some general arguments to this effect, such as the validity of the expansion improving as the 
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Figure 1: Comparison of asymptotic expressions (lower curve) versus sum over determinants (upper 
curve) for SU(5). 

dimensionality of the system increases. The only thing that can decide about this issue is a 
comparison with Monte Carlo simulations. In fact, it is known for some time that the Mean-Field 
predicts the order of the bulk phase transition and even the numerical value of the associated 
critical coupling correctly [1]. Note that this is a prediction already at the level of the MF 
background, i.e. without taking into account fluctuations. In [4] a more sophisticated comparison 
between the MF and Monte Carlo methods was performed on the anisotropic N = 2 lattice. The 
MF predicts, via the Wilson Loop (which can be obtained in the MF only by taking into account 
fluctuations), dimensional reduction via localization for 7 ~ 0.55 [2, 3]. This was verified by MC 
simulations that showed [4] that in the same regime of the anisotropy, Polyakov Loops fluctuate 
independently along four-dimensional hyperplanes. In [6] the mass of observables with 4d scalar 
and vector quantum numbers are computed with orbifold boundary conditions along the fifth 
dimension (these however should not interfere with the issue of the convergence of the MF 
expansion). Also these observables can be computed with the MF method only in the presence 
of fluctuations. The agreement was seen to be again good. All these results give us an increasing 
confidence that the MF is indeed a good description of the 5d system at the non-perturbative 
level. 

More specifically now, the computation of k$ from the gravity side suggested that the dual 
gauge theory, at leading order in a large N expansion a) must be not weakly coupled by the 
nature of the duality b) is in its 5d Coulomb phase since the static potential is of a 1/r 2 form 
c) sits in a regime of its phase diagram where cut-off effects are suppressed since in the result 
there is no sign of the presence of a gauge theory regulator. The first observation forced us to 
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0.06 1 1 1 1 1 1 1 1 1 

0.005 0.01 0.015 0.02 0.025 0.03 0.035 0.04 0.045 

1/L 



Figure 2: k 5 versus 1/L for SU{2) (lowest curve), 517(3), 5C/(4), SU(5), SU (10), 51/(15), 5[/(30) and 
517(100). 

move away from the 5d perturbative point where the coupling is expected to go to zero when 
the cut-off is removed. Moving into the interior of the phase diagram leads us eventually on the 
bulk phase transition because everywhere else physics is heavily cut-off dominated. The second 
observation on the other hand prohibits us from crossing the phase transition. The only choice 
we are left with is to be as near as possible to the phase transition, which is precisely where 
we have computed k$ from the gauge theory side. What is left to be addressed is why do we 
expect cut-off effects to be suppressed near the phase transition, or in lattice language, why do 
we generally expect the lattice spacing to decrease as ft — > (3 C . A milder version of this question, 
related to our quantity fcs is what is its lattice spacing dependence near the phase transition. It is 
not easy to answer definitively the former version of the question. What one can observe though 
from the MF computation [2] is that the scalar mass in units of lattice spacing amj decreases 
as the phase transition is approached. It is however a first order phase transition therefore most 
likely the lattice spacing remains finite and in fact one can see that ams can be pushed down 
to approximately 0.1 but not much further. This means that finite lattice spacing effects will 
be present even at the phase transition, even though not very large. Regarding k§, it has been 
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Figure 3: Intercept at L = oo versus 1/N along with a quadratic fit, predicting the value fcs = 0.0757 
for SU (oo). 



defined [8] by analogy to the non-abelian Coulomb's constant in four dimensions, where one can 
see that by construction it is a cut-off independent quantity. Here, our working assumption was 
that the lattice spacing dependence in k§ (via the product C2/3) cancels at least near the phase 
transition, just like in 4d and according to what the gravitational result suggests. The fact that 
as we move closer to the phase transition the value of k$ changes by very little, especially for 
large N, strongly supports this assumption. Notice that when L is increased while keeping q in 
Eq. (81) constant, takes us closer to the phase transition and thus changes the lattice spacing. 

We now discuss other possible sources for the observed 17% discrepancy. From the gravity 
side, in infinite four-dimensional volume, there may be a' and finite ./V corrections to the static 
potential. It would be interesting to see if it is possible to compute such corrections for the 
D4-brane background, which is the basis for the holographic computation of sect. 4. From 
the gauge theory side we also expect to have corrections. Already at leading order in the MF 
expansion, the computation of the coefficients A Z ,B Z involves truncations of series and of the 
ranks of the matrices the determinants of which determine the coefficients. Changing the order 
at which these are truncated would be a way to introduce error bars in the MF results for iV < 6. 
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Another would be to compute the coefficient C2 via local fits as in [3] as opposed to the global fits 
used here. This would require computing the Wilson Loops for increasing values of L (starting 
from approximately L = 200) and reading off each time the plateau value for C2- The variations 
in the plateau values would introduce additional error bars on the MF data. We did not do 
such an analysis here due to the lack of the necessary computing power. Finally, clearly there 
will be new effects if we change the order of the truncation of the MF expansion itself. Higher 
order effects will most certainly result in further corrections. It is hard to guess without further 
computations which of all possible corrections is mainly responsible for the discrepancy. 



6 Conclusion 



We developed the necessary formalism to perform the mean-field expansion to first non-trivial 
order, for five-dimensional, anisotropic SU (N) lattice gauge theories in a covariant gauge. We 
computed the mean- field background and then the propagator, the scalar observable and the 
Wilson Loop. The mass of the vector observable, being a geometric quantity, is the same as 
for SU(2). We then presented an application on the isotropic lattice. We computed Coulomb's 
constant in five dimensions and demonstrated that it becomes independent of N in the large N 
limit. At N = oo it converges to the value k$ = 0.0757, which amounts to approximately a 17% 
deviation from the value predicted by a holographic approach. 
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